Adiabatic Loading of Cold Bosons in Three-Dimensional Optical Lattices and 

Superfluid-Normal Phase Transition 
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We investigate the effects of the adiabatic loading of optical lattices to the temperature by applying 
the mean-field approximation to the three-dimensional Bose-Hubbard model at finite temperatures. 
We compute the lattice-height dependence of the isentropic curves for the given initial temperatures 
in case of the homogeneous system i.e., neglecting the trapping potential. Taking the unit of 
temperatures as the recoil energy, the adiabatic cooling/heating through superfiuid (SF) - normal 
(N) phase transition is clearly understood. It is found that the cooling occurs in SF phase while 
the heating occurs in N phase and the efficiency of adibatic cooling/heating is higher at higher 
temperatures. We also explain how its behavior can be understood from the lattice-hight dependence 
of dispersion relation in each phase. Furthermore, the connection of the adiabatic heating/cooling 
between the cases with/without the trapping potential is discussed. 



Recently ultracold atoms in optical lattices have been 
studied intensively both theoretically and experimen- 
tally (for the reviews, see 0, H U, 3)- Not only 
atomic, molecule, and optical (AMO) physics party, but 
also quantum information and condensed matter physics 
party come into this field and it has the possibilities of 
producing new kind of physics. From the viewpoint of 
quantum information, ultracold atoms in optical lattices 
can be used as one-way quantum computing Jp, 0, 0] 
and dynamical controlling of entanglement [8|, [9( under 
well-controlled conditions. Considering strongly corre- 
lated physics, this system offers the possibility of realizing 
various quantum lattice models like the Bose-Hubbard 
model, the Fermi-Hubbard model, and the Bose-Fermi 
Hubbard model, which have various rich quantum phases 
1- 

In order to investigate the above subjects, it is crucial 
to understand the lattice-height dependence of the tem- 
peratures. In the experiments, the temperature of Bosc 
gases is measured before inserting optical lattices. How- 
ever, the experimental method to investigate the tem- 
perature of Bose gases in optical lattices has not been es- 
tablished. Usually the loading process can be treated as 
adiabatic since the loading speed is very low. The behav- 
ior of the temperature of this system during the adiabatic 
loading of optical lattices is therefore of great interests, 
which has been studied in several papers 13, II , I2> IH • 



For the non-interacting Bose gases, the adiabatic cool- 
ing only occurs in the tight-binding regime (on the other 
hand, the adiabatic heating occurs when the thermal en- 
ergy lies in the first excited band). The mechanism of the 
adiabatic heating/cooling in this case can be understood 
in terms of the change of density of states [l0( • The adi- 
abatic loading including the interaction between atoms 
and the effect of the trapping potential has been studied 



for deep lattices. In that case, it is found that the adia- 
batic heating occurs due to increasing the Mott gap and 
the trapping effect induces the adiabatic compression and 
expansion which cause the adiabatic heating and cooling 
[TlL 12 . [l3j ]. On the other hand, the mechanism of the 
adiabatic cooling/heating through SF-N phase transition 
is not fully understood. 

In this paper, we examine the three-dimensional Bose- 
Hubbard model at finite temperatures within the mean- 
field approximation in order to investigate the effect of 
the adiabatic loading to the temperature of the system 
through SF phase to N phase. We show that the adi- 
abatic cooling occurs in SF phase while the adiabatic 
heating occurs in N phase. The number fluctuation con- 
serves during the adiabatic loading after crossing SF-N 
phase transition point, which yields that it is necessary 
to have ultracold temperatures at the phase transition 
point in order to obtain the system with very low num- 
ber fluctuation at deep optical lattices. We argue that 
the mechanism of the adiabatic cooling/heating is due 
to the dispersion relation in each phase. Finally, we will 
mention that in the case with the trapping potential, the 
mechanism of the adiabatic heating/cooling is essentially 
same as in the case without the trapping potential. 

The Hamiltonian for Bose atoms in optical lattices can 
be written as 



H = 
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V (x) 



d 3 x^{x)^{x)vb(x)i}(x) (1) 



where ip(x) is a field operator for Bose atoms and 
V (x) is the optical lattice potential. We consider three- 
dimensional optical lattices where V Q (x) has the form 

V {x) = V(s'm 2 kx + sin 2 ky + sin 2 kz). (2) 
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Here k = 2tt/X where A is the wave length of stand- 
ing wave laser forming optical lattices. The lattice con- 
stant is determined by a = A/2. The lattice height 
of optical lattices V is measured by the recoil energy 
Er = h 2 k 2 /2m where m is mass of the atom. We use the 
dimensionless lattice height s = V/Er. The binary in- 
teraction between atoms is approximated by s-wave scat- 
tering, which is charactrized by the scattering length a s . 

The Bosc- Hubbard Hamiltonian can be derived by ap- 
plying the tight-binding approximation to Eq. ([1]) [14l |. 
We expand the field operator by the Wannier function 

wq(x — Xi) as ffl(x) = J2i b\ Wo(x — Xi), where b\ is the 
destruction operator for a boson at a lattice site x\. We 
can rewrite Eq. fl} as 

h = -t J2 (p\ bj + h. c. ) + 1 J2 ^ fa - 1 ) - E ^ • ( 3 ) 

(ij) » i 

As usual, we consider the nearest-neighbor hopping and 
the on-site interaction. Note that we added the chemical 
potential /i in order to treat the grand-canonical ensem- 
ble. The relation between U, t and s under approximating 
the Wannier function as the Gaussian functioin can be 
found as = 5.97^s 88 and J- = 1.43s°- 98 e - 2 07 ^ 
[15| . In this paper, we consider 87 Rb and set m — 
1.44 x 1(T 25 [kg], a s = 545 [nm], and A = 852 [nm] 
as in the Greiner's experiment [l6j |. 

The validity of using the Bose-Hubbard model was ex- 
amined in Ref . 14( . The tight-binding approximation at 



finite temperatures implies that the energy scale includ- 
ing the thermal energy must be less than the gap energy 
between the lowest band and the second band. We cal- 
culated the thermal energy fc^T normalized by the gap 
energy A and confirmed fc^T/A <C 1 for the parameters 
we will use in this paper. 

Let us apply the mean-field approximation 

nam 

as 

+ bj) — (f> 2 where <fi = (b>) = (b) is taken 



bib, 



to be real. Then we can rewrite the Hamiltonian Eq.Q 
as sum of on-site Hamiltonians as H = Hi, and the 
Hamiltonian of the i-th site is 



Hi 



U 

—riiim - 1) 



zt(f>(b, 



(4) 



Here z denotes the coordinate number. In the actual con- 
ditions, we truncate the size of the Hilbert space of Hi 
by assuming the maximum number of particles which can 
be at one site is nt . We take large n t so that the trunca- 
tion effect on the calculated physical quantities is negleg- 
ible. Diagonalizing this Hamiltonian Hi under given U, 
t, and /i, we obtain the eigenstates and the correspond- 
ing eigenenergies. Then we can calculate the partition 
function Z and the Hclmholtz free energy F as functions 
of <j) for given temperatures. <j> is determined by the self- 
consistent equation dF/dcf) = 0. 

At first, we assume that the system is homogeneous 
i.e., neglecting the trapping potential. In order to inves- 
tigate the behavior of this system during the adiabatic 



loading of optical lattices, we calculate the entropy S/kB, 
the condensate density p s = <f> 2 , and the number fluctua- 
tion a = W (h 2 ) — (n) 2 under the fixed chemical potential 
ix as the occupation number p = 1 . Fig. [1] (a) shows the 
relation between the isentropic curves and the condensate 
density. Loading optical lattices adiabatically, the system 
goes along the isentropic curve which corresponds to the 
initial entropy. One finds that there exists two areas: the 
adiabatic cooling region and the adiabatic heating one. 
Former is the region where the system is cooled as the 
lattice height increases while the latter is that of where 
the system is heated as the lattice height increases. No- 
tice that the efficiency of the adiabatic cooling/heating is 
higher for higher initial entropy. Seeing the condensate 
density, one finds that SF phase, which is charcterized 
by p s > 0, exists in shallow lattices and at low temper- 
atures. The cooling region lies in SF phase, whereas the 
heating region lies in N phase where p s — 0. We also find 
that the efficiency of the adiabatic cooling/heating is not 
related with the amount of the condensate density. 

We plot the number fluctuation and the isentropic 
curves in Fig. Q] (b). For k B T/E R < 0.04, there exists 
a rectangle region with the very low number fluctuation, 
which is called as the thermal insulator in Ref. [l9| . In 
that region, the critical lattice height s c does not change 
much. It is necessary to satisfy S/ks < 0.01 in order 
to reach the thermal insulator. Note, however, that the 
number fluctuation is constant on the same isentropic 
curve in N phase. Thus, although the system is heated 
up, the number fluctuation does not change with loading 
optical lattices adiabatically after crossing SF-N phase 
transition point. 

One can understand the reason why the adiabatic cool- 
ing occurs in SF phase while the adiabatic heating does in 
N phase by investigating the dispersion relation in each 
phase (see Fig. [2]). The adiabatic loading means that 
the number of state is conserved during this process. As 
shown in Fig. [2] (a), the dispersion relation in SF phase 
exhibits the phonon-dispersion, where the sound veloc- 
ity decreases with increasing the lattice height from the 
initial state. Then the number of state which can be oc- 
cupied by the thermal energy increases if the temperature 
does not change. Thus the temperature must decrease in 
the final state in order to keep the number of state (i.e. 
the entropy) being constant. Fig.[2](b) shows that, in N 
phase, the energy gap increases as the lattice height in- 
creases from the initial state. Then the number of state 
which can be accesjsible by the thermal energy decreases 
if the temperature is same. Therefore the temperature 
rises to maintain the number of state being constant in 
the final state. In this way, the behavior of the adiabatic 
heating/cooling is determined by whether the dispersion 
relation increases or decreases from the initial state to 
the final state. 

Let us consider the system with the harmonic trap- 
ping potential Vt (r) = |w 2 r 2 , which usually exists in 
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FIG. 1: (Color online) (a) The isentropic curves and the con- 
densate density p s for the given lattice height s and temper- 
atures ksT/En. We choose the chemical potential p which 
satisfies p = 1. The number inside box on the isentropic 
curve indicates the value of the entropy S/ks- The system 
goes along the isentropic curve which corresponds to the ini- 
tial entropy. SF phase stands for the cooling region, whereas 
N phase stands for the heating region, (b) The isentropic 
curves added the density plot of the number fluctuation a. 
The isentropic curves is the same one as in (a). In SF phase, 
the number fluctuation is enhanced due to the quantum fluc- 
tuation. As loading optical lattices adiabatically, the num- 
ber fluctuation decreases and it takes the minimum value at 
the critical lattice height s c . The temperature increases after 
crosing s c whereas the number fluctuation is same. 
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FIG. 2: (Color online) Schematic picture of the dispersion re- 
lation and the thermal energy. Arrows indicate the direction 
of loading optical lattices. Its tail (head) corresponds to the 
initial (finial) state. (Dot: the thermal energy in the initial 
state (fcsT)ini, Short Dashed: the thermal energy in the fi- 
nal state (fcsT)fln, Long Dashed: The excitation spectrum in 
the initial state eini, Dot Dashed: The excitation spectrum 
in the final state efl n ) (a) SF phase. The sound velocity de- 
creases with increasing the lattice height. Thus the tempera- 
ture must decrease in order to keep the number of state, (b) 
N phase. The energy gap opens with loading optical lattices. 
Hence the temperature increase to maintain the entropy being 
a constant. 



experiments. Note that r is the distance from the center 
of the trapping potential. In the presence of the trap- 
ping potential, the system becomes inhomogenous and 
has the mixture of SF and N phases. In order to inves- 
tigate the behavior of the adiabatic heating/cooling in 
this case, we perform the same calculation as the homo- 
geneous case except inserting the local chemical poten- 
tial n(r) = no — Vt(r). We take the trapping frequency 



a- u> — 2ir x 24, the lattice size to be 65 3 , and the total num- 
2ber density N = 2 x 10 5 as in the Greiner's work [l6j |. 
Figs. (a) and (b) show qualitativley same results as 
homogeneous case. The adiabatic cooling occurs in the 
presence of condensate while the adiabatic heating oc- 
curs when there is no condensate. Therefore we can say 
that the behavior of adiabatic heating/cooling does not 
change essentially whether there is the trapping potential 
or not. We note, however, that the effects of the trap- 
ping potential is important quantitatively for realizing 
the strong correlated system as shown in Refs. 12, 
In Figs. [3] (c) and (d), we plot the spatial distributions 
of the occupation number, the condensate density, the 
entropy, and the number fluctuation along the adiabatic 
line Stotai/NkB — 0.3 indicated by two white dots in 
Figs. El (a) and (b). We see the high condensate density in 
Fig. [3] (c), while the wedding-cake structure is exhibited 
in Fig. [3] (d). Since we use local density approximation, 
further detailed studies are needed for the case with the 
trapping potential. 
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FIG. 3: (Color online) (a) The isentropic curves and the den- 
sity plot of the condensate density per particle, (b) The isoen- 
tropic curves and the density plot of the number fluctuation 
per particle, (c) and (d) Spatial distributions of p (circle, 
blue), p s (square, purple), S/ks (diamond, blown), and a 
(triangle, green) at (s,k B T/E R ) = (8.9, 0.085), (17.4, 0.045) 
respectively. This two (s, ksT / Er) are on the adiabatic line 
Stotai/Nks = 0.3, which are indicated by two white dots in 
(a) and (b). The condensate density exists in (c) while the 
wedding-cake structure is exbited in (d). 



In conclusion, we calculated the three-dimensional 
Bose-Hubbard model at finite temperatures within the 
mean-field approximation and investigated the effects of 
the adiabatic loading of optical lattices to the temper- 
ature. The lattice-height dependence of the isentropic 
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curves for given initial temperatures in case of the ho- 
mogeneous system was computed. We found that the 
cooling occurs in SF phase while the heating occurs in N 
phase and the efficiency of the adibatic cooling/heating 
is higher for higher temperatures. Its behavior is deter- 
mined by whether the dispersion relation of the system 
increases or decreases as loading optical lattices. Finally, 
we showed that the case with the trapping potential is 
essentially same as the homogeneous one. 

Note added. Recently we have become aware of a re- 
lated paper by Pollet et al. 20]. They studied the adi- 
abatic loading in one-dimensional and two-dimensional 
optical lattices by quantum monte carlo method. The 
effect of the trapping potential to the adiabatic loading 
was investigated with high accuracy in these lower di- 
mensional cases. 
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